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^ ■ ABSTRACT 

We measure the imprint of baryon acoustic oscillations (B AOs) in the galaxy clustering 
pattern at the highest redshift achieved to date, z — 0.6, using the distribution of A'' = 
132,509 emission-line galaxies in the WiggleZ Dark Energy Survey. We quantify BAOs 
using three statistics: the galaxy correlation function, power spectrum and the band- 
filtered estimator introduced by Xu et al. (2010). The results are mutually consistent, 
, corresponding to a 4.0% measurement of the cosmic distance-redshift relation at z = 

0.6 (in terms of the acoustic parameter "A(z)" introduced by Eisenstein et al. (2005) we 
find A{z — 0.6) — 0.452 ± 0.018). Both BAOs and power spectrum shape information 
contribute toward these constraints. The statistical significance of the detection of the 
acoustic peak in the correlation function, relative to a wiggle-free model, is 3.2-a. The 
ratios of our distance measurements to those obtained using BAOs in the distribution 
of Luminous Red Galaxies at redshifts z = 0.2 and z — 0.35 are consistent with a flat 
A Cold Dark Matter model that also provides a good fit to the pattern of observed 
fluctuations in the Cosmic Microwave Background (CMB) radiation. The addition 
of the current WiggleZ data results in a r:! 30% improvement in the measurement 
accuracy of a constant equation-of-state, w, using BAO data alone. Based solely on 
geometric BAO distance ratios, accelerating expansion {w < —1/3) is required with a 
probability of 99.8%, providing a consistency check of conclusions based on supernovae 
observations. Further improvements in cosmological constraints will result when the 
WiggleZ Survey dataset is complete. 
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1 INTRODUCTION 

The measurement of baryon acoustic oscillations (BAOs) in 
the large-scale clustering pattern of galaxies has rapidly be- 
come one of the most important observational pillars of the 
cosmological model. BAOs correspond to a preferred length 
scale imprinted in the distribution of photons and baryons 
by the propagation of sound waves in the rclativistic plasma 
of the early Universe (Peebles & Yu 1970, Sunyacv & Zel- 
dovitch 1970, Bond & Efstathiou 1984, Holtzman 1989, Hu 
& Sugiyama 1996, Eisenstein & Hu 1998). A full account 
of the early-universe physics is provided by Bashinsky &: 
Bertschinger (2001, 2002). In a simple intuitive description 
of the effect we can imagine an overdensity in the primor- 
dial dark matter distribution creating an overpressure in the 
tightly-coupled photon-baryon fluid and launching a spher- 
ical compression wave. At rcdshift z ~ 1000 there is a pre- 
cipitous decrease in sound speed due to recombination to a 
neutral gas and de-coupling of the photon-baryon fluid. The 
photons stream away and can be mapped as the Cosmic Mi- 
crowave Background (CMB) radiation; the spherical shell 
of compressed baryonic matter is frozen in place. The over- 
dense shell, together with the initial central perturbation, 
seeds the later formation of galaxies and imprints a preferred 
scale into the galaxy distribution equal to the sound hori- 
zon at the baryon drag epoch. Given that baryonic matter 
is secondary to cold dark matter in the clustering pattern, 
the amplitude of the effect is much smaller than the acoustic 
peak structure in the CMB. 

The measurement of BAOs in the pattern of late-time 
galaxy clustering provides a compelling validation of the 
standard picture that large-scale structure in today's Uni- 
verse arises through the gravitational amplification of per- 
turbations seeded at early times. The small amplitude of the 
imprint of BAOs in the galaxy distribution is a demonstra- 
tion that the bulk of matter consists of non-baryonic dark 
matter that does not couple to the relativistic plasma be- 
fore recombination. Furthermore, the preferred length scale 
- the sound horizon at the baryon drag epoch - may be pre- 
dicted very accurately by measurements of the CMB which 
yield the physical matter and baryon densities that control 
the sound speed, expansion rate and recombination time: 
the latest determination is 153.3 ± 2.0 Mpc (Komatsu et al. 
2009). Therefore the imprint of BAOs provide a standard 
cosmological ruler that can map out the cosmic expansion 
history and provide precise and robust constraints on the 
nature of the "dark energy" that is apparently dominating 
the current cosmic dynamics (Blake & Glazebrook 2003; Hu 
& Haiman 2003; Seo & Eisenstein 2003). In principle the 
standard ruler may be applied in both the tangential and 
radial directions of a gala:xy survey, yielding measures of 
the angular diameter distance and Hubble parameter as a 
function of rcdshift. 

The large scale and small amplitude of the BAOs im- 
printed in the galaxy distribution implies that galaxy red- 
shift surveys mapping cosmic volumes of order 1 Gpc^ with 
of order 10® galaxies are required to ensure a robust de- 
tection (Tegmark 1997, Blake & Glazebrook 2003, Blake et 
al. 2006). Gathering such a sample represents a formidable 
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observational challenge typically necessitating hundreds of 
nights of telescope time over several years. The leading such 
spectroscopic dataset in existence is the Sloan Digital Sky 
Survey (SDSS), which covers 8000 deg^ of sky containing a 
"main" r-band selected sample of lO'' galaxies with median 
redshift z « 0.1, and a Luminous Red Galaxy (LRG) exten- 
sion consisting of 10® galaxies but covering a significantly- 
greater cosmic volume with median redshift z w 0.35. Eisen- 
stein et al. (2005) reported a convincing BAO detection in 
the 2-point correlation function of the SDSS Third Data Re- 
lease (DR3) LRG sample aX z = 0.35, demonstrating that 
this standard-ruler measurement was self-consistent with the 
cosmological model established from CMB observations and 
yielding new, tighter constraints on cosmological parame- 
ters such as the spatial curvature. Pcrcival et al. (2010) 
undertook a power-spectrum analysis of the SDSS DR7 
dataset, considering both the main and LRG samples, and 
constrained the distance-redshift relation at both z — 0.2 
and z = 0.35 with ~ 3% accuracy in units of the standard 
ruler scale. Other studies of the SDSS LRG sample, pro- 
ducing broadly similar conclusions, have been performed by 
Huetsi (2006), Percival et al. (2007), Sanchez et al. (2009) 
and Kazin et al. (2010a). Some analyses have attempted 
to separate the tangential and radial BAO signatures in 
the LRG dataset, albeit with lower statistical significance 
(Gaztanaga et al. 2009, Kazin et al. 2010b). These studies 
built on earlier hints of BAOs reported by the 2-degree Field 
Galaxy Redshift Survey (Cole et al. 2005) and combinations 
of smaller datasets (Miller et al. 2001). 

This ambitious observational program to map out the 
cosmic expansion history with BAOs has prompted serious 
theoretical scrutiny of the accuracy with which we can model 
the BAO signature and the likely amplitude of systematic er- 
rors in the measurement. The pattern of clustering laid down 
in the high-redshift Universe is potentially subject to modu- 
lation by the non- linear scale-dependent growth of structure, 
by the distortions apparent when the signal is observed in 
redshift-space, and by the bias with which galaxies trace the 
underlying network of matter fluctuations. In this context 
the fact that the BAOs are imprinted on large, linear and 
quasi-linear scales of the clustering pattern implies that non- 
linear BAO distortions are relatively accessible to modelling 
via perturbation theory or numerical N-body simulations 
(Eisenstein, Seo & White 2007, Crocce & Scoccimarro 2008, 
Matsubara 2008). The leading-order effect is a "damping" of 
the sharpness of the acoustic feature due to the differential 
motion of pairs of tracers separated by 150 Mpc driven by 
bulk flows of matter. Effects due to galaxy formation and 
bias are confined to significantly smaller scales and are not 
expected to cause significant acoustic peak shifts. Although 
the non-linear damping of BAOs reduces to some extent the 
accuracy with which the standard ruler can be applied, the 
overall picture remains that BAOs provide a robust probe of 
the cosmological model free of serious systematic error. The 
principle challenge lies in executing the formidable galaxy 
redshift surveys needed to exploit the technique. 

In particular, the present ambition is to extend the 
relatively low-redshift BAO measurements provided by the 
SDSS dataset to the intermediate- and high-rcdsliift Uni- 
verse. Higher-redshift observations serve to further test the 
cosmological model over the full range of epochs for which 
dark energy apparently dominates the cosmic dynamics, can 



WiggleZ Survey: BAOs at z = 0.6 3 



probe greater cosmic volumes and therefore yield more ac- 
curate BAO measurements, and are less susceptible to the 
non-linear effects which damp the sharpness of the acoustic 
signature at low redshift and may induce low-amplitude sys- 
tematic errors. Currently, intermediate redshifts have only 
been probed by photometric-redshift surveys which have 
limited statistical precision (Blake et al. 2007, Padmanab- 
han et al. 2007). 

The WiggleZ Dark Energy Survey at the Australian As- 
tronomical Observatory (Drinkwater et al. 2010) was de- 
signed to provide the next-generation spectroscopic BAO 
dataset following the SDSS, extending the distance-scale 
measurements across the intermediate-redshift range up to 
z = 0.9 with a precision of mapping the acoustic scale com- 
parable to the SDSS LRG sample. The survey, which began 
in August 2006, completed observations in January 2011 
and has obtained of order 200,000 redshifts for UV-bright 
emission-line galaxies covering of order 1000 deg^ of equato- 
rial sky. Analysis of the full dataset is ongoing. In this paper 
we report intermediate results for a subset of the WiggleZ 
sample with effective redshift z = 0.6. 

BAOs are a signature present in the 2-point clustering 
of galaxies. In this paper we analyze this signature using a 
variety of techniques: the 2-point correlation function, the 
power spectrum, and the band-filtered estimator recently 
proposed by Xu et al. (2010) which amounts to a band- 
filtered correlation function. Quantifying the BAO measure- 
ment using this range of techniques increases the robustness 
of our results and gives us a sense of the amplitude of sys- 
tematic errors induced by our current methodologies. Us- 
ing each of these techniques we measure the angle-averaged 
clustering statistic, making no attempt to separate the tan- 
gential and radial components of the signal. Therefore we 
measure the "dilation scale" distance Dv{z) introduced by 
Eisenstein et al. (2005) which consists of two parts physi- 
cal angular-diameter distance, Da{z), and one part radial 
proper-distance, cz/H{z): 



Dv{z) = 



1/3 



(1) 



This distance measure reflects the relative importance of 
the tangential and radial modes in the angle-averaged BAO 
measurement (Padmanabhan & White 2008), and reduces 
to proper distance in the low-redshift limit. Given that a 
measurement of Dv{z) is correlated with the physical mat- 
ter density Qmh^ which controls the standard ruler scale, we 
extract other distilled parameters which are far less signifi- 
cantly correlated with Qmh^, namely: the acoustic parame- 
ter A{z) as introduced by Eisenstein et al. (2005); the ratio 
dz = rs{zd) / Dv{z), which quantifies the distance scale in 
units of the sound horizon at the baryon drag epoch, rs(z;d); 
and 1/Rz which is the ratio between Dy (z) and the distance 
to the CMB last-scattering surface. 

The structure of this paper is as follows. The Wig- 
gleZ data sample is introduced in Section 2, and we then 
present our measurements of the galaxy correlation function, 
power spectrum and band-filtered correlation function in 
Sections 3, 4 and 5, respectively. The results of these differ- 
ent methodologies are compared in Section 6. In Section 7 we 
state our measurements of the BAO distance scale at z = 0.6 
using various distilled parameters, and combine our result 
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Figure 1. The probability distribution of galaxy redsliifts in each 
of the WiggleZ regions used in our clustering analysis, together 
with the combined distribution. Differences between individual 
regions result from variations in the galaxy colour selection crite- 
ria depending on the available optical imaging (Drinkwater et al. 
2010). 



with other cosmological datasets in Section 8. Throughout 
this paper we assume a fiducial cosmological model which 
is a flat ACDM Universe with matter density parameter 
Sim ~ 0.27, baryon fraction Qh/^m = 0.166, Hubble pa- 
rameter h = 0.71, primordial index of scalar perturbations 
Jig = 0.96 and redshift-zero normalization erg — 0.8. This 
fiducial model is used for some of the intermediate steps 
in our analysis but our final cosmological constraints are, 
to first-order at least, independent of the choice of fiducial 
model. 



2 DATA 

The WiggleZ Dark Energy Survey at the Anglo Australian 
Telescope (Drinkwater et al. 2010) is a large-scale galaxy 
redshift survey of bright emission-line galaxies mapping a 
cosmic volume of order 1 Gpc^ over the redshift interval 
2 < 1. The survey has obtained of order 200,000 redshifts 
for UV-selected galaxies covering of order 1000 deg^ of equa- 
torial sky. In this paper we analyze the subset of the WiggleZ 
sample assembled up to the end of the lOA semester (May 
2010). We include data from six survey regions in the red- 
shift range 0.3 < z < 0.9 - the 9-hr, 11-hr, 15-hr, 22-hr, 1-hr 
and 3-hr regions - which together constitute a total sample 
of A'' = 132,509 galaxies. The redshift probability distribu- 
tions of the galaxies in each region are shown in Figure 1. 

The selection function for each survey region was deter- 
mined using the methods described by Blake et al. (2010) 
which model effects due to the survey boundaries, incom- 
pleteness in the parent UV and optical catalogues, incom- 
pleteness in the spectroscopic follow-up, systematic varia- 
tions in the spectroscopic redshift completeness across the 
AAOmega spectrograph, and variations of the galaxy red- 
shift distribution with angular position. The modelling pro- 
cess produces a series of Monte Carlo random realizations of 
the angle/redshift catalogue in each region, which are used 
in the correlation function estimation. By stacking together 
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a very large number of these random realizations we deduced 
the 3D window function grid used for power spectrum esti- 
mation. 



3 CORRELATION FUNCTION 
3.1 Measurements 

The 2-point correlation function is a common method for 
quantifying the clustering of a population of galaxies, in 
which the distribution of pair separations in the dataset 
is compared to that within random, unclustered catalogues 
possessing the same selection function (Peebles 1980). In 
the context of measuring baryon acoustic oscillations, the 
correlation function has the advantage that the expected 
signal of a preferred clustering scale is confined to a single, 
narrow range of separations around 105 Mpc. Further- 
more, small-scale non-linear effects, such as the distribution 
of galaxies within dark matter haloes, do not influence the 
correlation function on these large scales. One disadvantage 
of this statistic is that measurements of the large-scale cor- 
relation function are prone to systematic error because they 
are very sensitive to the unknown mean density of the galaxy 
population. However, such "integral constraint" effects re- 
sult in a roughly constant offset in the large-scale correla- 
tion function, which does not introduce a preferred scale 
that could mimic the BAG signature. 

In order to estimate the correlation function of each 
WiggleZ survey region we first placed the angle/redshift cat- 
alogues for the data and random sets on a grid of co-moving 
co-ordinates, assuming a flat ACDM model with matter den- 
sity fim = 0.27. We then measured the redshift-space 2-point 
correlation function ^(s) for each region using the Landy- 
Szalay (1993) estimator: 

DD{s) - DR(s) + RR{s) 

^^'^ = Rm ' 

where DD{s), DR{s) and RR{s) are the data-data, data- 
random and random-random weighted pair counts in sep- 
aration bin s, each random catalogue containing the same 
number of galaxies as the real dataset. In the construction 
of the pair counts each data or random galaxy i is assigned 
a weight Wi — 1/(1 + riiPo), where n; is the survey num- 
ber density [in Mpc~^] at the location of the ith galaxy, 
and Po = 5000 /i"^ Mpc^ is a characteristic power spectrum 
amplitude at the scales of interest. The survey number den- 
sity distribution is established by averaging over a large en- 
semble of random catalogues. The DR and RR pair counts 
are determined by averaging over 10 random catalogues. We 
measured the correlation function in 20 separation bins of 
width 10 Mpc between 10 and 180 /i"^ Mpc, and de- 
termined the covariance matrix of this measurement using 
lognormal survey realizations as described below. We com- 
bined the correlation function measurements in each bin for 
the different survey regions using inverse- variance weighting 
of each measurement (we note that this procedure produces 
an almost identical result to combining the individual pair 
counts). 

The combined correlation function is plotted in Figure 
2 and shows clear evidence for the baryon acoustic peak at 
separation ~ 105 Mpc. The effective redshift ZcS of the 
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Figure 2. The combined redshift-space correlation function ^(s) 
for WiggleZ survey regions, plotted in tlie combination ^(s) 
where s is the co-moving redshift-space separation. The best- 
fitting clustering model (varying Qm'i^, o and b^) is overplot- 
ted as the solid line. We also show as the dashed line the corre- 
sponding "no- wiggles" reference model, constructed from a power 
spectrum with the same clustering amplitude but lacking baryon 
acoustic oscillations. 

correlation function measurement is the weighted mean red- 
shift of the galaxy pairs entering the calculation, where the 
redshift of a pair is simply the average {zi + Z2)/2, and the 
weighting is wi'W2 where Wi is defined above. We determined 
Zcff for the bin 100 < s < 110 Mpc, although it does not 
vary significantly with separation. For the combined Wig- 
gleZ survey measurement, we found z^g = 0.60. 

We note that the correlation function measurements are 
corrected for the effect of redshift blunders in the WiggleZ 
data catalogue. These are fully quantified in Section 3.2 of 
Blake et al. (2010), and can be well-approximated by a scale- 
independent boost to the correlation function amplitude of 
(1 — fb)^^, where fb ~ 0.05 is the redshift blunder fraction 
(which is separately measured for each WiggleZ region). 

3.2 Uncertainties : lognormal realizations and 
covariance matrix 

We determined the covariance matrix of the correlation func- 
tion measurement in each survey region using a large set of 
lognormal realizations. Jack-knife errors, implemented by di- 
viding the survey volume into many sub-regions, are a poor 
approximation for the error in the large-scale correlation 
function because the pair separations of interest are usually 
comparable to the size of the sub-regions, which are then 
not strictly independent. Furthermore, because the WiggleZ 
dataset is not volume-limited and the galaxy number den- 
sity varies with position, it is impossible to define a set of 
sub-regions which are strictly equivalent. 

Lognormal realizations are relatively cheap to generate 
and provide a reasonably accurate galaxy clustering model 
for the linear and quasi-linear scales which are important for 
the modelling of baryon oscillations (Coles & Jones 1991). 
We generated a set of realizations for each survey region 
using the method described in Blake & Glazebrook (2003) 
and Glazebrook & Blake (2005). In brief, we started with 
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Figure 3. The amplitude of tiie cross-correlation dj / y^CuCjj 
of the covariance matrix dj for the correlation function mea- 
surement plotted in Figure 2, determined using lognormal real- 
izations. 



a model galaxy power spectrum Pmod(fc) consistent with 
the survey measurement. We then constructed Gaussian re- 
alizations of overdensity fields Sair) sampled from a sec- 
ond power spectrum Pcik) « Pmod(fc) (defined below), in 
which real and imaginary Fourier amplitudes are drawn 
from a Gaussian distribut ion with zero mean and stan- 
dard deviation \/ PG{k)/2. A lognormal overdensity field 
5LN(fO ~ exp((5G) — 1 is then created, and is used to pro- 
duce a galaxy density field Pg(r) consistent with the survey 
window function W{r): 



Pg{r)(xW{r) [l + 5LN(r)], 



(3) 



where the constant of proportionality is fixed by the size 
of the final dataset. The galaxy catalogue is then Poisson- 
sampled in cells from the density field Pg(r). We note that 
the input power spectrum for the Gaussian overdensity 
field, Pcik), is constructed to ensure that the final power 
spectrum of the lognormal overdensity field is consistent 
with Pmod(fc). This is achieved using the relation between 
the correlation functions of Gaussian and lognormal fields, 
^G(r-) = ln [l + Cmod«]. 

We determined the covariance matrix between bins i 
and j using the correlation function measurements from a 
large ensemble of lognormal realizations: 



(4) 



where the angled brackets indicate an average over the re- 
alizations. Figure 3 displays the final covariance matrix re- 
sulting from combining the different WiggleZ survey regions 
in the form of a correlation matrix dj / CuCjj. The mag- 
nitude of the first and second off-diagonal elements of the 
correlation matrix is typically 0.6 and 0.4, respectively. We 
find that the jack-knife errors on scales of 100 Mpc typ- 
ically exceed the lognormal errors by a factor of « 50%, 
which we can attribute to an over-estimation of the number 
of independent jack-knife regions. 



3.3 Fitting the correlation function : template 
model and simulations 

In this Section we discuss the construction of the template 
fiducial correlation function model ^fld, galaxy (s) which we fit- 
ted to the WiggleZ measurement. When fitting the model we 
vary a scale distortion parameter a, a linear normalization 
factor and the matter density Q,uih? which controls both 
the overall shape of the correlation function and the stan- 
dard ruler sound horizon scale. Hence we fitted the model 



^mod(s) = b ^fld,galaxy(Q! S). 



(5) 



The probability distribution of the scale distortion param- 
eter a, after marginalizing over fl,mh^ and b^, gives the 
probability distribution of the distance variable Dv{zcb) ~ 
ct Dv,fid{zoH) where z^s — 0.6 for our sample (Eisenstein 
et al. 2005, Padmanabhan & White 2008). Dy, defined by 
Equation 1, is a composite of the physical angular-diameter 
distance Da{z) and Hubble parameter H{z) which govern 
tangential and radial galaxy separations, respectively, where 
I>v,fld(zcff) = 2085.4 Mpc. 

We note that the measured value of Dv resulting from 
this fitting process will be independent (to first order) of 
tfic fiducial cosmological model adopted for the conversion 
of galaxy redshifts and angular positions to co-moving co- 
ordinates. A change in Dvfid would result in a shift in the 
measured position of the acoustic peak. This shift would 
be compensated for by a corresponding offset in the best- 
fitting value of a, leaving the measurement of Dv = a Dvfid 
unchanged (to first order). 

An angle-averaged power spectrum P{k) may be con- 
verted into an angle-averaged correlation function ^(s) using 
the spherical Hankel transform 



C(s) = ^ / dkk'Pik) 



sin (ks) 



ks 



(6) 



In order to determine the shape of the model power spec- 
trum for a given flmh^, we first generated a linear power 
spectrum Phik) using the fitting formula of Eisenstein & 
Hu (1998). This yields a result in good agreement with a 
GAME linear power spectrum (Lewis, Challinor & Lasenby 
2000), and also produces a wiggle- free reference spectrum 
Prcf(k) which possesses the same shape as Phik) but with 
the baryon oscillation component deleted. This reference 
spectrum is useful for assessing the statistical significance 
with which we have detected the acoustic peak. We fixed 
the values of the other cosmological parameters using our 
fiducial model h = 0.71, fib^^ = 0.0226, = 0.96 and 
(78 = 0.8. Our choices for these parameters are consistent 
with the latest fits to the Cosmic Microwave Background 
radiation (Komatsu et al. 2009). 

We then corrected the power spectrum for quasi-linear 
effects. There are two main aspects to the model: a damping 
of the acoustic peak caused by the displacement of matter 
due to bulk fiows, and a distortion in the overall shape of 
the clustering pattern due to the scale-dependent growth of 
structure (Eisenstein, Seo & White 2007, Grocce & Scocci- 
marro 2008, Matsubara 2008). We constructed our model in 
a similar manner to Eisenstein et al. (2005). We first incorpo- 
rated the acoustic peak smoothing by multiplying the power 
spectrum by a Gaussian damping term g{k) = exp {—k^al): 



apcd{k) = g{k) PL(k) + [1 - g{k)] P,^t(k), 



(7) 
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where the inclusion of the second term maintains the same 
small-scale clustering amplitude. The magnitude of the 
damping can be modelled using perturbation theory (Crocce 
& Scoccimarro 2008) as 



67r2 



PL{k)dk, 



(8) 



where / — f2m(z)'''^^ is the growth rate of structure. In our 
fiducial cosmological model, fl^h^ = 0.136f, we find = 
4.5 Mpc. We checked that this value was consistent with 
the allowed range when a-v was varied as a free parameter 
and fitted to the data. 

Next, we incorporated the non- linear boost to the clus- 
tering power using the fitting formula of Smith et al. (2003). 
However, we calculated the non- linear enhancement of power 
using the input no-wiggles reference spectrum rather than 
the full linear model including baryon oscillations: 



damped, NL 



(fc) = 



Prcf(fc) 



X -Pdampcd (^) • 



(9) 



Equation 9 is then transformed into a correlation function 
Cdamped,NL(s) usiug Equatiou 6. 

The final component of our model is a scale-dependent 
galaxy bias term B{s) relating the galaxy correlation func- 
tion appearing in Equation 5 to the non-linear matter cor- 
relation function: 



Cfid, galaxy (s) = -B(s) ^damped, NL (s) , 



(10) 



where we note that an overall constant normalization has 
already been separated in Equation 5 so that B{s) — >■ 1 at 
large s. 

We determined the form of B{s) using halo catalogues 
extracted from the GiggleZ dark matter simulation. This 
TV-body simulation has been generated specifically in sup- 
port of WiggleZ survey science, and consists of 2160^ par- 
ticles evolved in a 1 h~'^ Gpc^ box using a WMAP5 cos- 
mology (Komatsu et al. 2009). We deduced B{s) using the 
non-linear redshift-space halo correlation functions and non- 
linear dark-matter correlation function of the simulation. 
We found that a satisfactory fitting formula for the scale- 
dependent bias over the scales of interest is 



(11) 



We performed this procedure for several contiguous sub- 
sets of 250,000 halos rank-ordered by their maximum cir- 
cular velocity (a robust proxy for halo mass). The best- 
fitting parameters of Equation 11 for the subset which best 
matches the large-scale WiggleZ clustering amplitude are 
So — 0.32 Mpc, 7 = —1.36. We note that the magnitude 
of the scale-dependent correction from this term is ~ 1% 
for a scale s ~ 10 Mpc, which is far smaller than the 
~ 10% magnitude of such effects for more strongly-biased 
galaxy samples such as Luminous Red Galaxies (Eisenstein 
et al. 2005). This greatly reduces the potential for system- 
atic error due to a failure to model correctly scale-dependent 
galaxy bias effects. 



3.4 Extraction of Dy 

We fitted the galaxy correlation function template model 
described above to the WiggleZ survey measurement, vary- 
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Figure 4. The scale-dependent correction to the non-linear real- 
space dark matter correlation function for haloes with maximum 
circular velocity Vmax ~ 125 km s~^, which possess the same 
amplitude of large-scale clustering as WiggleZ galaxies. The green 
line is the ratio of the real-space halo correlation function to the 
real-space non-linear dark matter correlation function. The red 
line is the ratio of the redshift-space halo correlation function 
to the real-space halo correlation function. The black line, the 
product of the red and green lines, is the scale-dependent bias 
correction B(s) which we fitted with the model of Equation 11, 
shown as the dashed black line. The blue line is the ratio of the 
real-space non-linear to linear correlation function. 



ing the matter density Qmh^, the scale distortion parame- 
ter a and the galaxy bias fe^. Our default fitting range was 
10 < s < 180 /i"^ Mpc (following Eisenstein et al. 2005), 
where 10 is an estimate of the minimum scale of va- 
lidity for the quasi-linear theory described in Section 3.3. 
In the following, we assess the sensitivity of the parameter 
constraints to the fitting range. 

We minimized the statistic using the full data co- 
variance matrix, assuming that the probability of a model 
was proportional to exp(— x^/2)- The best-fitting parame- 
ters were flmh'^ = 0.132 ± 0.011, a = 1.075 ± 0.055 and 
= 1.21±0.11, where the errors in each parameter are pro- 
duced by marginalizing over the remaining two parameters. 
The minimum value of is 14.9 for 14 degrees of freedom 
(17 bins minus 3 fitted parameters), indicating an accept- 
able fit to the data. In Figure 2 we compare the best-fitting 
correlation function model to the WiggleZ data points. The 
results of the parameter fits are summarized for ease of ref- 
erence in Table 1. 

Our measurement of the scale distortion parameter a 
may be translated into a constraint on the distance scale 
Dv = aDvfid = 2234.9 ± 115.2 Mpc, corresponding to a 
5.2% measurement of the distance scale at z = 0.60. This 
accuracy is comparable to that reported by Eisenstein et al. 
(2005) for the analysis of the SDSS DR3 LRG sample at z ^ 
0.35. Figure 5 compares our measurement of the distance- 
redshift relation with those from the LRG samples analyzed 
by Eisenstein et al. (2005) and Percival et al. (2010). 

The 2D probability contours for the parameters Qmh^ 
and Dv{z = 0.6), marginalizing over b^, are displayed in Fig- 
ure 6. Following Eisenstein et al. (2005) we indicate three 
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Figure 5. Measurements of the distance-redshift relation using 
the BAO standard ruler from LRG samples (Eisenstein et al. 2005, 
Percival et al. 2010) and the current WiggleZ analysis. The results 
are compared to a fiducial flat ACDM cosmological model with 
matter density Qrn = 0.27. 

degeneracy directions in this parameter space. The first di- 
rection (the dashed line in Figure 6) corresponds to a con- 
stant measured acoustic peak separation, i.e. rs{zd) / Dv{z = 
0.6) = constant. We used the fitting formula quoted in Per- 
cival et al. (2010) to determine rs(zd) as a function of Quih^ 
(given our fiducial value of Qbh^ = 0.0226); we find that 
ra{zd) = 152.6 Mpc for our fiducial cosmological model. The 
second degeneracy direction (the dotted line in Figure 6) cor- 
responds to a constant measured shape of a Cold Dark Mat- 
ter power spectrum, i.e. Dv{z — 0.6) x Q^nh^ ~ constant. In 
such models the matter transfer function at recombination 
can be expressed as a function of g = k/Q^nh^ (Bardeen et 
al. 1986). Given that changing Dv corresponds to a scal- 
ing of oc Dv.fid/ Dv , we recover that the measured power 
spectrum shape depends on Dvi^mh^- The principle de- 
generacy axis of our measurement lies between these two 
curves, suggesting that both the correlation function shape 
and acoustic peak information are driving our measurement 
of Dv- The third degeneracy direction we plot (the dash- 
dotted line in Figure 6), which matches our measurement, 
corresponds to a consta nt value of the acoustic parameter 
A{z) = Dv {z)\/ QuiHq jcz introduced by Eisenstein et al. 
(2005). We present our fits for this parameter in Section 7. 

In Figure 6 we also show probability contours resulting 
from fits to a restricted range of separations s > 30 and 
50 hr^ Mpc. In both cases the contours become significantly 
more extended and the long ajcis shifts into alignment with 
the case of the acoustic peak alone driving the fits. The 
restricted fitting range no longer enables us to perform an 
accurate determination of the value of ilmft^ from the shape 
of the clustering pattern alone. 

3.5 Significance of the acoustic peak detection 

In order to assess the importance of the baryon acoustic 
peak in constraining this model, we repeated the parameter 
fit replacing the model correlation function with one gener- 
ated using a "no- wiggles" reference power spectrum Prof (A:), 
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Figure 6. Probability contours of the physical matter density 
f2m/i^ and distance scale Dv(z = 0.6) obtained by fitting to the 
WiggleZ survey combined correlation function ^(s). Results are 
compared for different ranges of fitted scales Sj„in < s < 180 
Mpc. The (black solid, red dashed, blue dot-dashed) contours cor- 
respond to fitting for Sniin — (10, 30, 50)/i-i Mpc, respectively 
The heavy dashed and dotted lines are the degeneracy directions 
which are expected to result from fits involving respectively just 
the acoustic oscillations, and just the shape of a pure CDM power 
spectrum. The heavy dash-dotted line represents a constant value 
of the acoustic "yl" parameter introduced by Eisenstein et al. 
(2005), which is the parameter best-measured by our correlation 
function data. The solid circle represents the location of our fidu- 
cial cosmological model. The two contour levels in each case en- 
close regions containing 68% and 95% of the likelihood. 

which possesses the same amplitude and overall shape as 
the original matter power spectrum but lacks the baryon 
oscillation features (i.e., we replaced Phiji) with Pref(fc) in 
Equation 7) . The minimum value obtained for the x'^ statis- 
tic in this case was 25.0, indicating that the model containing 
baryon oscillations was favoured by Ax^ = 10.1. This corre- 
sponds to a detection of the acoustic peak with a statistical 
significance of 3.2-(j. Furthermore, the value and error ob- 
tained for the scale distortion parameter in the no-wiggles 
model was ol = 0.80 ± 0.17, representing a degradation of 
the error in a by a factor of three. This also suggests that 
the acoustic peak is important for establishing the distance 
constraints from our measurement. 

As an alternative approach for assessing the significance 
of the acoustic peak, we changed the fiducial baryon density 
to Sib = and repeated the parameter fit. The minimum 
value obtained for the statistic was now 22.7 and the 
value and marginalized error determined for the scale dis- 
tortion parameter was a = 0.80 ± 0.12, re-affirming the sig- 
nificance of our detection of the baryon wiggles. 

If we restrict the correlation function fits to the range 
50 < s < 130 Mpc, further reducing the influence of the 
overall shape of the clustering pattern on the goodness-of- 
fit, we find that our fiducial model has a minimum x^ ~ 5.9 
(for 5 degrees of freedom) and the "no-wiggles" reference 
spectrum produces a minimum x'^ = 13.1. Even for this 
restricted range of scales, the model containing baryon os- 
cillations was therefore favoured by Ax^ ~ 7.2. 
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3.6 Sensitivity to the clustering model 

In this Section we investigate tlie systematic dependence 
of our measurement of Dv{z = 0.6) on the model used to 
describe the quasi-linear correlation function. We considered 
five modelling approaches proposed in the literature: 

• Model 1: Our fiducial model described in Section 3.3 
following Eisenstein et al. (2005), in which the quasi-linear 
damping of the acoustic peak was modelled by an expo- 
nential factor g{k) — exp (— fc^a^), is determined from 
linear theory via Equation 8, and the small-scale power was 
restored by adding a term [l — g{k)] multiplied by the wiggle- 
free reference spectrum (Equation 7). 

• Model 2: No quasi-linear damping of the acoustic peak 
was applied, i.e. (t„ = 0. 

• Model 3: The term restoring the small-scale power, [1 — 
(7(fe)]Prof (fc) in Equation 7, was omitted. 

• Model 4: fdamped(fc) in Equation 7 was generated us- 
ing Equation 14 of Eisenstein, Seo & White (2006), which 
implements different damping coefficients in the tangential 
and radial directions. 

• Model 5: The quasi-linear matter correlation function 
was generated using Equation 10 of Sanchez et al. (2009), 
following Crocce & Scoccimarro (2008), which includes the 
additional contribution of a "mode-coupling" term. We set 
the coefficient Amc = 1 in this equation (rather than intro- 
duce an additional free parameter). 

Figure 7 compares the measurements of Dv{z = 0.6) 
from the correlation function data, marginalized over Qmh^ 
and b^, assuming each of these models. The agreement 
amongst the best-fitting measurements is excellent, and the 
minimum statistics imply a good fit to the data in each 
case. We conclude that systematic errors associated with 
modelling the correlation function are not significantly af- 
fecting our results. The error in the distance measurement 
is determined by the amount of damping of the acoustic 
peak, which controls the precision with which the standard 
ruler may be applied. The lowest distance error is produced 
by Model 2 which neglects damping; the greatest distance 
error is associated with Model 4, in which the damping is 
enhanced along the line-of-sight (see Equation 13 in Eisen- 
stein, Seo & White 2006). 



4 POWER SPECTRUM 

4.1 Measurements and covariance matrix 

The power spectrum is a second commonly-used method 
for quantifying the galaxy clustering pattern, which is com- 
plementary to the correlation function. It is calculated us- 
ing a Fourier decomposition of the density field in which 
(contrary to the correlation function) the maximal signal- 
to-noise is achieved on large, linear or quasi-linear scales 
(at low wavenumbers) and the measurement of small-scale 
power (at high wavenumbers) is limited by shot noise. How- 
ever, also in contrast to the correlation function, small-scale 
effects such as shot noise infiuence the measured power at 
all wavenumbers, and the baryon oscillation signature ap- 
pears as a series of decaying harmonic peaks and troughs at 
different wavenumbers. In aesthetic terms this diffusion of 
the baryon oscillation signal is disadvantageous. 
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Figure 7. Measurements of Dy{z = 0.6) from the galaxy corre- 
lation function, marginalizing over flmh^ and , comparing five 
different models for the quasi-linear correlation function as de- 
tailed in the text. The measurements are consistent, suggesting 
that systematic modelling errors are not significantly affecting our 
results. 

We estimated the galaxy power spectrum for each sep- 
arate WiggleZ survey region using the direct Fourier meth- 
ods introduced by Feldman, Kaiser & Peacock (1994; FKP). 
Our methodology is fully described in Section 3.1 of Blake 
et al. (2010); we give a brief summary here. Firstly we map 
the angle-redshift survey cone into a cuboid of co-moving 
co-ordinates using a fiducial fiat ACDM cosmological model 
with matter density Q^n ~ 0.27. We gridded the catalogue in 
cells using nearest grid point assignment ensuring that the 
Nyquist frequencies in each direction were much higher than 
the Fourier wavenumbers of interest (we corrected the power 
spectrum measurement for the small bias introduced by this 
gridding). We then applied a Fast Fourier transform to the 
grid, optimally weighting each pixel by 1/(1 + nPo), where 
n is the galaxy number density in the pixel (determined us- 
ing the selection function) and Po = 5000 Mpc^ is a 
characteristic power spectrum amplitude. The Fast Fourier 
transform of the selection function is then used to construct 
the final power spectrum estimator using Equation 13 in 
Blake et al. (2010). The measurement is corrected for the 
effect of redshift blunders using Monte Carlo survey simula- 
tions as described in Section 3.2 of Blake et al. (2010). We 
measured each power spectrum in wavenumber bins of width 
0.01 h Mpc"^ between k = and 0.3 h Mpc"\ and deter- 
mined the covariance matrix of the measurement in these 
bins by implementing the sums in Fourier space described 
by FKP (see Blake et al. 2010 equations 20-22). The FKP 
errors agree with those obtained from lognormal realizations 
within 10% at all scales. 

In order to detect and fit for the baryon oscillation sig- 
nature in the WiggleZ galaxy power spectrum, we need to 
stack together the measurements in the individual survey re- 
gions and redshift slices. This requires care because each sub- 
region possesses a different selection function, and therefore 
each power spectrum measurement corresponds to a differ- 
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Table 1. Results of fitting a three-parameter model {fl^h^,a,b^) to WiggleZ measurements of four different clustering statistics for 
various ranges of scales. The top four entries, above the horizontal line, correspond to our fiducial choices of fitting range for each statistic. 
The fitted scales a are converted into measurements of Dy and two BAO distilled parameters, A and rs{z^)/Dy, which are introduced 
in Section 7. The final column lists the measured value of Dy when the parameter f2mfe^ is left fixed at its fiducial value and only the bias 

is marginalized. We recommend using A{z = 0.6) as measured by the correlation function ^(s) for the scale range 10 < s < 180 
Mpc, highlighted in bold, as the most appropriate WiggleZ measurement for deriving BAO constraints on cosmological parameters. 



Statistic 


Scale range 




Dviz = 0.6) 


A{z = 0.6) 


rs(z,,)/Dv{z = 0.6) 


Dv-(z = 0.6) 








[Mix-] 






fixing O^Ji^ 




10 < s < ISOh-^ Mpc 


0.1,32 ±0.011 


2234.9 ± 115.2 


0.452 ±0.018 


0.0692 ± 0.0033 


2216.5 ± 78.9 


P{k) [full] 


0.02 < < 0.2/1 Mpc-i 


0.134 ±0.008 


2160.7 ± 132.3 


0.440 ± 0.020 


0.0711 ±0.0038 


2141.0 ± 97.5 


P(fc) [wiggles] 


0.02 <k <0.2h Mpc-l 


0.163 ±0.017 


2135.4 ± 156.7 


0.461 ±0.030 


0.0699 ±0.0045 


2197.2 ± 119.1 


wo (r) 


10 < r < 180/1^1 Mpc 


0.130 ±0.011 


2279.2 ± 142.4 


0.456 ±0.021 


0.0680 ±0.0037 


2238.2 ± 104.6 



5(s) 30 < s < 180/j-i Mpc 0.166 ±0.014 2127.7 ± 127.9 0.475 ± 0.025 0.0689 ± 0.0031 2246.8 ± 102.6 

5(s) 50 < s < 180/j-i Mpc 0.164 ±0.016 2129.2 ± 140.8 0.474 ± 0.025 0.0690 ± 0.0031 2240.1 ± 104.7 

P(A;) [full] 0.02 < A; < O.l/i Mpc-i 0.150 ±0.020 2044.7 ±253.0 0.441 ± 0.034 0.0733 ± 0.0073 2218.1 ± 128.4 

P(fe) [full] 0.02 < fc < 0.3/1 Mpc" 1 0.137±0.007 2132.1 ± 109.2 0.441 ± 0.017 0.0716 ± 0.0033 2148.9 ± 79.9 

P(*;) [wiggles] 0.02 < fc < 0.1 /i Mpc" ^ 0.160 ±0.020 2240.7 ±235.8 0.466 ± 0.034 0.0678 ± 0.0070 2277.9 ± 187.5 

P(/c) [wiggles] 0.02 < fe < 0.3/1 Mpc" 1 0.161 ±0.019 2114.5 ± 132.4 0.455 ± 0.026 0.0706 ± 0.0037 2171.4 ±98.0 

wo{r) 30 < r < 180/i-i Mpc 0.127 ±0.018 2288.8 ± 157.3 0.455 ± 0.027 0.0681 ± 0.00,37 2251.6 ± 111.7 

wo{r) 50 < r < 180/1-1 Mpc 0.164 ±0.016 2190.0 ± 146.2 0.466 ± 0.023 0.0673 ± 0.0036 2282.1 ± 109.8 



ent convolution of the underlying power spectrum model. 
Furthermore the non-linear component of the underlying 
model varies with redshift, due to non-Iincar evolution of 
the density and velocity power spectra. Hence the observed 
power spectrum in general has a systematically-different 
slope in each sub-region, which implies that the baryon os- 
cillation peaks lie at slightly different wavenumbers. If we 
stacked together the raw measurements, there would be a 
significant washing-out of the acoustic peak structure. 

Therefore, before combining the measurements, we 
made a correction to the shape of the various power spectra 
to bring them into alignment. We wish to avoid spuriously 
enhancing the oscillatory features when making this correc- 
tion. Our starting point is therefore a fiducial power spec- 
trum model generated from the Eisenstein & Hu (1998) "no- 
wiggles" reference linear power spectrum, which defines the 
fiducial slope to which we correct each measurement. Firstly, 
we modified this reference function into a redshift-space non- 
linear power spectrum, using an empirical redshift-space dis- 
tortion model fitted to the two-dimensional power spectrum 
split into tangential and radial bins (see Blake et al. 2011a). 
The redshift-space distortion is modelled by a coherent-flow 
parameter /? and a pairwise velocity dispersion parameter 
Ov, which were fitted independently in each of the redshift 
slices. We convolved this redshift-space non-linear reference 
power spectrum with the selection function in each sub- 
region, and our correction factor for the measured power 
spectrum is then the ratio of this convolved function to the 
original real-space linear reference power spectrum. After 
applying this correction to the data and covariance ma^ 
trix wc combined the resulting power spectra using inverse- 
variance weighting. 

Figures 8 and 9 respectively display the combined power 
spectrum data, and that data divided through by the com- 
bined no-wiggles reference spectrum in order to reveal any 
signature of acoustic oscillations more clearly. We note that 
there is a significant enhancement of power at the position of 



the first harmonic, k « 0.075 h Mpc~^. The other harmonics 
are not clearly detected with the current dataset, although 
the model is nevertheless a good statistical fit. Figure 10 
displays the final power spectrum covariance matrix, result- 
ing from combining the different WiggleZ survey regions, 
in the form of a correlation matrix Cjj/y^ CuCjj. We note 
that there is very little correlation between separate 0.01 h 
Mpc^i power spectrum bins. 

We note that our method for combining power spec- 
trum measurements in different sub-regions only corrects for 
the convolution effect of the window function on the overall 
power spectrum shape, and does not undo the smoothing 
of the BAO signature in each window. We therefore expect 
the resulting BAO detection in the combined power spec- 
trum may have somewhat lower significance than that in 
the combined correlation function. 



4.2 Extraction of Dy 

We investigated two separate methods for fitting the scale 
distortion parameter to the power spectrum data. Our first 
approach used the whole shape of the power spectrum in- 
cluding any baryonic signature. We generated a template 
model non-linear power spectrum Psdik) parameterized by 
Qmh^ , which we took as Equation 9 in Section 3.3, and fitted 
the model 

Pmod(fc) = b^Pfid{k/a), (12) 

where a now appears in the denominator (as opposed to 
the numerator of Equation 5) due to the switch from real 
space to Fourier space. As in the case of the correlation func- 
tion, the probability distribution of a, after marginalizing 
over Qmh and , can be connected to the measurement of 
Dv{zes)- We determined the effective redshift of the power 
spectrum estimate by weighting each pixel in the selection 
function by its contribution to the power spectrum error: 
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Figure 8. The power spectrum obtained by staclcing measure- 
ments in different WiggleZ survey regions using the method de- 
scribed in Section 4.1. The best-fitting power spectrum model 
(varying 0,^^? , 01 and b^) is overplotted as the solid line. We 
also show the corresponding "no-wiggles" reference model as the 
dashed line, constructed from a power spectrum with the same 
clustering amplitude but lacking baryon acoustic oscillations. 
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Reference 
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Figure 9. The combined WiggleZ survey power spectrum of Fig- 
ure 8 divided by the smooth reference spectrum to reveal the 
signature of baryon oscillations more clearly. We detect the first 
harmonic peak in Fourier space. 



ZcB 



(13) 



where ng{x) is the galaxy number density in each grid cell 
X and Pg is the characteristic galaxy power spectrum am- 
plitude, which we evaluated at a scale k = 0.1 h Mpc~^. We 
obtained an effective redshift z^h = 0.583. In order to enable 
comparison with the correlation function fits we applied the 
best-fitting value of a at 2: = 0.6. 

Our second approach to fitting the power spectrum 
measurement used only the information contained in the 
baryon oscillations. We divided the combined WiggleZ 



Figure 10. The amplitude of the cross-correlation Cij / ^Ca Cjj 
of the covariance matrix dj for the power spectrum measure- 
ment, determined using the FKP estimator. The amplitude of 
the off-diagonal elements of the covariance matrix is very low. 



power spectrum data by the corresponding combined no- 
wiggles reference spectrum, and when fitting models we di- 
vided each trial power spectrum by its corresponding refer- 
ence spectrum prior to evaluating the statistic. 

We restricted our fits to Fourier wavescales 0.02 < k < 
0.2 h Mpc~^, where the upper limit is an estimate of the 
range of reliability of the quasi-linear power spectrum mod- 
elling. We investigate below the sensitivity of the best-fitting 
parameters to the fitting range. For the first method, fitting 
to the full power spectrum shape, the best-fitting parame- 
ters and 68% confidence ranges were fim/i^ = 0.134 ±0.008 
and a = 1.050 ± 0.064, where the errors in each parameter 
are produced by marginalizing over the remaining two pa- 
rameters. The minimum value of was 12.4 for 15 degrees 
of freedom (18 bins minus 3 fitted parameters), indicating 
an acceptable fit to the data. We can convert the constraint 
on the scale distortion parameter into a measured distance 
Dv{z = 0.6) = 2160.7 ± 132.3 Mpc. The 2D probability dis- 
tribution of Qmh^ and Dv{z — 0.6), marginalizing over b^, 
is displayed as the solid contours in Figure 11. In this Figure 
we reproduce the same degeneracy lines discussed in Section 
3.4, which are expected to result from fits involving just the 
acoustic oscillations and just the shape of a pure CDM power 
spectrum. We note that the long axis of our probability con- 
tours is oriented close to the latter line, indicating that the 
acoustic peak is not exerting a strong influence on flts to the 
full WiggleZ power spectrum shape. Comparison of Figure 
11 with Figure 6 shows that flts to the WiggleZ galaxy cor- 
relation function are currently more influenced by the BAOs 
than the power spectrum. This is attributable to the signal 
being stacked at a single scale in the correlation function, in 
this case of a moderate BAO detection. 

For the second method, fltting to just the baryon os- 
cillations, the best-fltting parameters and 68% confldence 
ranges were Clmh'^ = 0.163 ± 0.017 and a = 1.000 ± 0.073. 
Inspection of the 2D probability contours of Q^h^ and a, 
which are shown as the dotted contours in Figure 11, indi- 
cates that a signiflcant degeneracy has opened up parallel 
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Figure 11. Probability contours of the physical matter den- 
sity Oni/i^ and distance scale Dy{z = 0.6) obtained by fitting 
to the WiggleZ survey combined power spectrum. Results are 
compared for different ranges of fitted scales 0.02 < k < femax 
and methods. The (red dashed, black solid, green dot-dashed) 
contours correspond to fits of the full power spectrum model 
for fcniax = (0.1, 0.2, 0.3) /i Mpc-l, respectively. The blue dot- 
ted contours result from fitting to the power spectrum divided 
by a smooth no- wiggles reference spectrum (with fcmax = 0.2 h 
Mpc~^). Degeneracy directions and likelihood contour levels are 
plotted as in Figure 6. 



to the line of constant apparent BAO scale (as expected). 
Increasing f2m/i^ decreases the standard ruler scale, but the 
positions of the acoustic peaks may be brought back into 
line with the data by applying a lower scale distortion pa- 
rameter a. Low values of Qmh^ are ruled out because the 
resulting amplitude of baryon oscillations is too high (given 
that Q]^h^ is fixed). We also plot in Figure 11 the probabil- 
ity contours resulting from fitting different ranges of Fourier 
scales fc < 0.1 ft Mpc"^ and k < 0.3/i Mpc'^ The 68% con- 
fidence regions generated for these different cases overlap. 

We assessed the significance with which acoustic fea- 
tures are detected in the power spectrum using a method 
similar to our treatment of the correlation function in Sec- 
tion 3. We repeated the parameter fit for (fim/i^, a, b^) using 
the "no-wiggles" reference power spectrum in place of the 
full model power spectrum. The minimum value obtained 
for the statistic in this case was 15.8, indicating that the 
model containing baryon oscillations was favoured by only 
Ax^ = 3.3. This is consistent with the direction of the long 
axis of the probability contours in Figure 11, which suggests 
that the baryon oscillations are not driving the fits to the 
full power spectrum shape. 



5 BAND-FILTERED CORRELATION 
FUNCTION 

5.1 Measurements and covariance matrix 

Xu et al. (2010) introduced a new statistic for the measure- 
ment of the acoustic peak in configuration space, which they 
describe as an advantageous approach for band-filtering the 
information. They proposed estimating the quantity 



(14) 



where ^(s) is the 2-point correlation function as a function 
of separation s and 



otherwise 



0<x <1 (15) 
(16) 



in terms of a; = (s/r)^. This filter localizes the acoustic in- 
formation in a single feature at the acoustic scale in a similar 
manner to the correlation function, which is beneficial for se- 
curing a robust, model-independent detection. However, the 
form of the filter function W{x) advantageously reduces the 
sensitivity to small-scale power (which is difficult to model 
due to non-linear effects) and large-scale power (which is 
difficult to measure because it is subject to uncertainties 
regarding the mean density of the sample), combining the 
respective advantages of the correlation function and power 
spectrum approaches. 

Xu et al. (2010) proposed that wo{r) should be esti- 
mated as a weighted sum over galaxy pairs i 



i— 1 



(17) 



where Si is the separation of pair i, \ii is the cosine of the 
angle of the separation vector to the line-of-sight, Nd is the 
number of data galaxies, and no is the average galaxy den- 
sity which we simply define as Nd /V where V is the volume 
of a cuboid enclosing the survey cone. The function (j>{s, jj.) 
describes the edge effects due to the survey boundaries and 
is normalized so that the number of random pairs in a bin 
(s — >■ s -|- ds, n ^ fi + d/i) is 



RR{s, ^) = 2-KnDNDS (l>{s, ^) ds dfi, 



(18) 



where (p = 1 for a uniform, infinite survey. We determined 
the function (j){s, fi) used in Equation 17 by binning the pair 
counts RR{s,fi) for many random sets in fine bins of s and 
/i. We then fitted a parameterized model 



(19) 



to the result in bins of s, and used the coefficients a„(s) to 
generate the value of (f) for each galaxy pair. 

We note that our Equation 17 contains an extra factor 
of 2 compared to Equation 12 of Xu et al. (2010) because we 
define the quantity as a sum over unique pairs, rather than 
all pairs. We also propose to modify the estimator to intro- 
duce a "DR" term by analogy with the correlation function 
estimator of Equation 2, in order to correct for the distribu- 
tion of data galaxies with respect to the boundaries of the 
sample: 



wo{r) = DDmtcicd(r) - DRmtcrcdir), 



(20) 



where DRmtcrcd{r) is estimated using Equation 17, but sum- 
ming over data-random pairs and excluding the initial factor 
of 2. We used our lognormal realizations to determine that 
this modified estimator of Equation 20 produces a result 
with lower bias and variance compared to Equation 17. We 
used Equation 20 to measure the band filtered correlation 
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Best-fitting model 
Reference 




50 100 150 

r / In"' Mpc 

Figure 12. The band-filtered correlation function wo{r) for the 
combined WiggleZ survey regions, plotted in the combination 
r-^wo{r). The best-fitting clustering model (varying Qmh^, a and 
b^) is overplotted as the solid line. We also show the corresponding 
"no-wiggles" reference model, constructed from a power spectrum 
with the same clustering amplitude but lacking baryon acoustic 
oscillations. We note that the high covariance of the data points 
for this estimator implies that (despite appearances) the solid line 
is a good statistical fit to the data. 



function of each WiggleZ region for 17 values of r spaced by 
10 Mpc between 15 and 175 /i"^ Mpc. 

We determined the covariance matrix dj of our esti- 
mator using the ensemble of lognormal realizations for each 
survey region. We note that for our dataset the amplitude of 
the diagonal errors \/Cii determined by lognormal realiza- 
tions is typically ~ 5 times greater than jack-knife errors and 
~ 3 times higher than obtained by evaluating Equation 13 
in Xu et al. (2010) which estimates the covariance matrix in 
the Gaussian limit. Given the likely drawbacks of jack-knife 
errors (the lack of independence of the jack-knife regions on 
large scales) and Gaussian errors (which fail to incorporate 
the survey selection function), the lognormal errors should 
provide by far the best estimate of the covariance matrix for 
this measurement. 

We constructed the final measurement of the band- 
filtered correlation function by stacking the individual mea- 
surements in different survey regions with inverse-variance 
weighting. Figure 12 displays our measurement. We detect 
clear evidence of the expected dip in wo(r) at the acoustic 
scale. Figure 13 displays the final covariance matrix of the 
band-filtered correlation function resulting from combining 
the different WiggleZ survey regions in the form of a corre- 
lation matrix Cij / y/CuCjj . We note that the nature of the 
wo{r) estimator, which depends on the correlation function 
at all scales s < r, implies that the data points in different 
bins of r are highly correlated, and the correlation coefficient 
increases with r. At the acoustic scale, neighbouring 10 
Mpc bins are correlated at the ~ 85% level and bins spaced 
by 20 h-^ Mpc are correlated at a level of ~ 55%. 



5.2 Extraction of Dy 

We determined the acoustic scale from the band-filtered 
correlation function by constructing a template function 
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Figure 13. The amplitude of the cross-correlation dj / y^Cu Cjj 
of the covariance matrix dj for the band-filtered correlation func- 
tion measurement plotted in Figure 12, determined using lognor- 
mal realizations. 



Wo. fid. galaxy (j') in the Same style as Section 3.3 and then fit- 
ting the model 



lUO,mod('') = 6^ WO, fid, galaxy (a r). 



(21) 



We determined the function wo, Ad, galaxy (f) by applying the 
transformation of Equation 14 to the template galaxy cor- 
relation function ^fld, galaxy (s) defined in Section 3.3, as a 
function of Qmh^ ■ 

The best-fitting parameters to the band-filtered corre- 
lation function are D.n^h'^ = 0.130 ±0.011, a = 1.100 ±0.069 
and — 1.32 ± 0.13, where the errors in each parameter 
are produced by marginalizing over the remaining two pa- 
rameters. The minimum value of is 10.5 for 14 degrees of 
freedom (17 bins minus 3 fitted parameters), indicating an 
acceptable fit to the data. Our measurement of the distor- 
tion parameter may be translated into a constraint on the 
distance scale Dv{z = 0.6) = aDyfid = 2279.2 ±142.4 Mpc, 
corresponding to a 6.2% measurement of the distance scale 
at z = 0.6. Probability contours of Dv{z — 0.6) and Qmh^ 
are overplotted in Figure 14. 

In Figure 12 we compare the best fitting band-filtered 
correlation function model to the WiggleZ data points (not- 
ing that the strong covariance between the data gives the 
mis-leading impression of a poor fit). We overplot a second 
model which corresponds to our best-fit parameters but for 
which the no- wiggle reference power spectrum has been used 
in place of the full power spectrum. If we fit this no-wiggles 
model varying Qmh^ , a and 6^ we find a minimum value 
of = 19.0, implying that the model containing acoustic 
features is favoured by Ax^ = 8.5. 



6 COMPARISON OF CLUSTERING 
STATISTICS 

The distance-scale measurements of Dv{z = 0.6) using the 
four different clustering statistics applied in this paper are 
compared in Figure 14. All four statistics give broadly con- 
sistent results for the measurement of Dy and Slm/i^, with 
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Figure 14. A comparison of tiie probability contours of Q^h^ 
and Dy{z = 0.6) resulting from fitting different clustering statis- 
tics measured from the WiggleZ survey: the correlation func- 
tion for scale range 10 < s < 180 ft,- 1 Mpc (the solid black 
contours), the band-filtered correlation function for scale range 
10 < r < 180 Mpc (the dashed red contours), the full power 
spectrum shape for scale range 0.02 < fc < 0.2h Mpc~^ (the dot- 
dashed green contours), and the power spectrum divided by a 
"no- wiggles" reference spectrum (the dotted blue contours). De- 
generacy directions and likelihood contour levels are plotted as in 
Figure 6. 



For each of the clustering statistics determined above, the 
measurement of Dv is significantly correlated with the mat- 
ter density Q.mh^ which controls both the shape of the clus- 
tering pattern and the length-scale of the standard ruler (see 
Figure 14). It is therefore useful to re-cast these BAO mea- 
surements in a manner less correlated with Q.uih^ and more 
representative of the observable combination of parameters 
constrained by the BAOs. These "distilled parameters" are 
introduced and measured in this Section. 



7.1 CMB information 

The length-scale of the BAO standard ruler and shape of 
the linear clustering pattern is calibrated by Cosmic Mi- 
crowave Background data. The cosmological information 
contained in the CMB may be conveniently encapsulated 
by the Wilkinson Microwave Anisotropy Probe (WMAP) 
"distance priors" (Komatsu et al. 2009). We use the 7-year 
WMAP results quoted in Komatsu et al. (2011). 

Firstly, the CMB accurately measures the characteris- 
tic angular scale of the acoustic peaks 6a = rs(z*)/(l + 
Zt)DA{z*), where rs{zt) is the size of the sound horizon at 
last scattering and Da{z*) is the physical angular-diameter 
distance to the decoupling surface. This quantity is conven- 
tionally expressed as a characteristic acoustic index: 



Ia = tv/Oa 



n{l + z,)^^4^ = 302.09 ± 0.76. 



(22) 



significant overlap between the respective 68% confidence 
regions in this parameter space. This agreement suggests 
that systematic measurement errors in these statistics are 
not currently dominating the WiggleZ BAO fits. 

There are some differences in detail between the results 
derived from the four statistics. Fits to the galaxy power 
spectrum are currently dominated by the power spectrum 
shape rather than the BAOs, such that the degeneracy direc- 
tion lies along the line of constant apparent turnover scale. 
Fitting to only the "wiggles" in Fourier space gives weaker 
constraints on the distance scale which (unsurprisingly) lie 
along the line of constant apparent BAO scale. 

The correlation function and band-filtered correlation 
function yield very similar results (with the constraints from 
the standard correlation function being slightly stronger). 
Their degeneracy direction in the {Q,mh^,Dv) parameter 
space lies between the two degeneracy directions previously 
mentioned, implying that both the BAO scale and correla- 
tion function shape are infiuencing the result. The slightly 
weaker constraint on the distance scale provided by the 
band-filtered correlation function compared to the standard 
correlation function is likely due to the suppression of infor- 
mation on small and large scales by the compensated filter, 
which is designed to reduce potential systematic errors in 
modelling the shape of the clustering pattern. 

For our cosmological parameter fits in the remainder of 
the paper we used the standard correlation function as our 
default choice of statistic. The correlation function provides 
the tightest measurement of the distance scale from our cur- 
rent dataset and encodes the most significant detection of 
the BAO signal. 



The complete CMB likelihood is well-reproduced by com- 
bining this measurement of £a with the "shift parameter" 
defined by 



7^; 



(1 + z^)Da{z,) = 1.725 ± 0.018 



(23) 



and the redshift of recombination (using the fitting function 
given as equations 66-68 in Komatsu et al. 2009) 



1091.3 ±0.91. 



(24) 



The inverse covariance matrix for {(.AjTI, z*) is given as Ta- 
ble 10 in Komatsu et al. (2011) and is included in our cos- 
mological parameter fit. 



7.2 Measuring A{z) 

As noted in Eisenstein et al. (2005) and discussed in Section 
3.4 above, the parameter combination 



Aiz) 



Dv{z)y^, 



ff2 



(25) 



which we refer to as the "acoustic parameter", is particu- 
larly well-constrained by distance fits which utilize a com- 
bination of acoustic oscillation and clustering shape infor- 
mation, since in this situation the degeneracy direction of 
constant A{z) lies approximately perpendicular to the mi- 
nor axis of the measured (Dy,f2m/i^) probability contours. 
Conveniently, A{z) is also independent of Ho (given that 
Dv oc l/Ho). Figure 15 displays the measurements result- 
ing from fitting the parameter set {A,Qjnh^ ,b^) to the four 
WiggleZ clustering statistics and marginalizing over . The 
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Figure 15. A comparison of the results of fitting different Wig- 
gleZ clustering measurements in the same style as Figure 14, ex- 
cept that we now fit for the parameter A{z = 0.6) (defined by 
Equation 25) rather than Dv{z = 0.6). 



results of the parameter fits are displayed in Table 1; the cor- 
relation function yields A{z = 0.6) = 0.452 ±0.018 (i-e. with 
a measurement precision of 4.0%). For this clustering statis- 
tic in particular, the correlation between measurements of 
A(z) and ^l^h? is very low. Given that the CMB provides a 
very accurate determination of Q.mh? (via the distance pri- 
ors) we do not use the WiggleZ determination of Q,uih^ in 
our cosmological parameter fits, but just use the marginal- 
ized measurement of A[z). 

The measurement of A(z) involves the assumption of 
a model for the shape of the power spectrum, which we 
parameterize by ilmft^- Essentially the full power spectrum 
shape, rather than just the BAOs, is being used as a stan- 
dard ruler, although the two features combine in such a way 
that A and Q,uih^ are uncorrelated. However, given that a 
model for the full power spectrum is being employed, we re- 
fer to these results as large-scale structure ("LSS") rather 
than BAO constraints, where appropriate. 



7.3 Measuring dz 

In the case of a measurement of the BAOs in which the 
shape of the clustering pattern is marginalized over, the 
{Dv,^mh?) probability contours would lie along a line of 
constant apparent BAO scale. Hence the extracted distances 
are measured in units of the standard-ruler scale, which 
may be conveniently quoted using the distilled parameter 
dz = ra{zd)/Dv{z) where rs{zd) is the co-moving sound 
horizon size at the baryon drag epoch. In contrast to the 
acoustic parameter, dz provides a purely geometric distance 
measurement that does not depend on knowledge of the 
power spectrum shape. The information required to com- 
pare the observations to theoretical predictions also varies 
between these first two distilled parameters: the prediction 
of dz requires prior information about h (or Quih^), whereas 
the prediction of A{z) does not. We also fitted the parameter 
set (do. 6, ^imh^, b^) to the four WiggleZ clustering statistics. 
The results of the parameter fits are displayed in Table 1; 



the correlation function yields do.e ~ 0.0692 ± 0.0033 (i.e. 
with a measurement precision of 4.8%). We note that, given 
the WiggleZ fits are in part driven by the shape of the power 
spectrum as well as the BAOs, there is a weak residual cor- 
relation between do. 6 and flmh^- 

When calculating the theoretical prediction for this pa- 
rameter we obtained the value of Vs {zd) for each cosmological 
model tested using Equation 6 of Eisenstein & Hu (1998), 
which is a fitting formula for rs(zd) in terms of the values 
of Qmh^ and Q^^li^ . In our analysis we fixed Q]^h^ — 0.0226 
which is consistent with the measured CMB value (Komatsu 
et al. 2009); we find that marginalizing over the uncertainty 
in this value does not change the results of our cosmological 
analysis. 

We note that rs{zd) is determined from the matter and 
baryon densities in units of Mpc (not Mpc), and thus a 
fiducial value of h must also be used when determining dz 
from data (we chose h — 0.71). However, the quoted observa- 
tional result do. 6 ~ 0.0692 ± 0.0033 is actually independent 
of h. Adoption of a different value of h would result in a 
shifted standard ruler scale [in units of h"^ Mpc] and hence 
shifted best-fitting values of a and Dy in such a way that 
dz is unchanged. However, although the observed value of 
dz is independent of the fiducial value of h, the model fitted 
to the data still depends on h as remarked above. 



7.4 Measuring Rz 

The measurement of dz may be equivalently expressed as a 
ratio of the low-redshift distance Dv{z) to the distance to 
the last-scattering surface, exploiting the accurate measure- 
ment of £a provided by the CMB. We note that the value of 
Rz depends on the behaviour of dark energy between red- 
shift z and recombination, whereas a constraint derived from 
dz only depends on the properties of dark energy at redshifts 
lower than z. Taking the product of dz and Ia/tt approx- 
imately cancels out the dependence on the sound horizon 
scale: 



1/Rz = Ia dz/Ti 



{l + z.) 



Da[z,) rsjzd) 
rs{z,) Dv(z) 



(26) 



- (l + ..):^x 1.044. 
Dv{z) 

The value 1.044 is the ratio between the sound horizon at 
last scattering and at the baryon drag epoch. Although this 
is a model-dependent quantity, the change in redshift be- 
tween recombination and the end of the drag epoch is driven 
by the relative number density of photons and baryons, 
which is a feature that does not change much across the 
range of viable cosmological models. Combining our mea- 
surement of do.6 = 0.0692 ± 0.0033 from the WiggleZ corre- 
lation function fit with £a ~ 302.09 ± 0.76 (Komatsu et al. 
2011) we obtain l/Ro.a = 6.65 ± 0.32. 



7.5 Measuring distance ratios 

Finally, we can avoid the need to combine the BAO fits with 
CMB measurements by considering distance ratios between 
the different redshifts at which BAO detections have been 
performed. Measurements of Dv{z) alone are dependent on 
the fiducial cosmological model and assumed standard-ruler 
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scale: an efficient way to measure Dviz2) /Dv{zi) is by cal- 
culating dzi/dz2 , which is independent of the value of rs{zd). 
Percival et al. (2010) reported BAO fits to the SDSS LRG 
sample in two correlated redshift bins do. 2 = 0.1905±0.0061 
and do. 35 — 0.1097 ± 0.0036 (with correlation coefficient 
0.337). Ratioing do. 2 with the independent measurement of 
do.6 = 0.0692 ± 0.0033 from the WiggleZ correlation func- 
tion and combining the errors in quadrature we find that 
Dv{0-6)/Dv{0.2) = 2.753 ± 0.158. Percival et al. (2010) 
report Dv{0.35)/Dv{0-2) = 1.737 ± 0.065 (where in this 
latter case the error is slightly tighter than obtained by 
adding errors in quadrature because of the correlation be- 
tween Dv{0.2) and Dv{0.35)). These two distance ratio 
measurements are also correlated by the common presence 
of Dv^(0.2) in the denominators; the correlation coefficient 
is 0.313. 

7.6 Comparison of distilled parameters 

Figure 16 compares the measurements of the distilled pa- 
rameters introduced in this Section, using both WiggleZ and 
LRG data, to a series of cosmological models varying either 
Qui or the dark energy equation-of-state, w, relative to a 
fiducial model with Q^n = 0.27 and w = —1. These plots 
have all been normalized to the fiducial model and plotted 
on the same scale so that the level of information contained 
in the different distilled parameters can be compared. For 
presentational purposes in Figure 16 we converted the Per- 
cival et al. (2010) measurements of dz into constraints on 
the acoustic parameter A{z), using their fiducial values of 
rs{zd) = 154.7 Mpc and Om/i^ = 0.1296 and employing the 
same fractional error in the distilled parameter. 

We note that A{z) and the ratio of Dv{z) values are 
the only distilled parameters that are independent of CMB 
data: the measurement of 1/Rz uses £a from the CMB, and 
whilst the observed value of dz is independent of the CMB, 
we need to marginalize over Q,rnh^ or h in order to compare 
it with theoretical models. Figure 17 illustrates the differing 
information encapsulated by each of these distilled parame- 
ters by plotting likelihood contours of {Qm, Qa) for w = —1 
and (ilm, w) for ilk = fitted to these data. The significant 
differences in the resulting likelihood contours exemplify the 
fact that in some cases we have not used all the information 
in the galaxy data (e.g. Dv{z) / Dv{0.2) uses only the ratio of 
BAO scales, neglecting information from the absolute scale 
of the standard ruler), while in other cases we have already 
included significant information from the CMB (e.g. 1/Rz). 
This illustrates that we must be careful to include additional 
CMB data in a self-consistent manner, not double-counting 
the information. This plot also serves to illustrate the mild 
tension between the BAO and CMB results, which can be 
seen from the fact that the 1/Rz constraints are offset from 
the centre of the A{z) constraints. The WiggleZ and LRG 
distilled BAO parameters used in these fits are summarized 
for convenience in Table 2. 



8 COSMOLOGICAL PARAMETER 
MEASUREMENTS 

In this Section we present cosmological parameter fits to 
the distilled BAO parameters measured above. We consider 
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Figure 16. Measurements of the distilled BAO parameters ex- 
tracted from WiggleZ and LRG datasets as a function of red- 
shift. From upper to lower we plot Dy (in Mpc), the acous- 
tic parameter ^4(2), = rs{zd)/ Dy{z), 1/Rz = Ia^z/'^ and 
Dv{z)/ Dv{z = 0.2). The left-hand column shows a set of differ- 
ent cosmological models varying Om with [Q^, w] = [0, —1] in each 
case. The right-hand column displays a range of models varying 
w with [Cm, = [0-27, 0] in each case. 

two versions of the standard cosmological model. The first 
version is the standard ACDM model in which dark energy 
is a cosmological constant with equation-of-state w = —1 
and spatial curvature is a free parameter; we fit for Q.rn and 
the cosmological constant density f^A. The second version is 
the fiat uiCDM model in which spatial curvature is fixed at 
f2k = but the equation of state of dark energy is a free 
parameter; we fit for fim and w. 

Unless otherwise stated, we fitted these cosmological 
models to the WiggleZ measurement of A{z = 0.6) combined 
with the measurements of dz from LRG samples at 2; = 0.2 
and z = 0.35 (Percival et al. 2010). The distilled param- 
eter A{z) (measured from the galaxy correlation function) 
is the most appropriate choice for quantifying the WiggleZ 
BAO measurement because it is uncorrelated with S7mft'^, as 
demonstrated by Figure 15. The parameter dz provides the 
best representation of the Percival et al. (2010) BAO data. 
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Figure 17. Likelihood contours (l-cr and 2-a) derived from model fits to the different distilled parameters which may be used to 
encapsulate the BAO results from LRG and WiggleZ data. In each case the contours represent the combination of the three redshift 
bins for which we have BAO data, z = [0.2,0.35,0.6]. The thick solid lines with grey shading show the A{z) parameter constraints, 
which are the most appropriate representation of WiggleZ data. The three sets of green lines show constraints from dz assuming three 
different priors: green solid lines include a prior fim/*^ = 0.1326 it 0.0063 (Komatsu et al. 2009, as used by Percival et al. 2010); green 
dotted lines marginalize over a flat prior of 0.5 < h < 1.0; and green dashed lines marginalize over a Gaussian prior of h = 0.72 ± 0.03. 
Red dashed lines show the 1/Rz constraints, whilst fits to the ratios of BAO distances Dy{z)/ Dy{0.2) are shown by blue dash-dot 
lines. The left-hand panel shows results for curved cosmological-constant universes parameterized by (Omi ^^a) ^^<^ the right-hand panel 
displays results for fiat dark-energy universes parameterized by (fim, w). Comparisons between these contours reveal the differing levels of 
information encoded in each distilled parameter. By combining each type of distilled parameter with the CMB data in a correct manner, 
self-consistent results should be achieved. 



Table 2. Measurements of the distilled BAO parameters at redshifts z = 0.2, 0.35 and 0.6 from LRG and WiggleZ data, which are used 
in our cosmological parameter constraints. The LRG measurements at z = 0.2 and z = 0.35 are correlated with coefficient 0.337 (Percival 
et al. 2010). The two distance ratios rfo.2/<io.35 E'nd rfo.2/c'o.6 correlated with coefficient 0.313. The different measurements of 1/Rz 
are correlated by the common presence of the £a variable and by the covariance between do. 2 and do. 35- Our default cosmological fits 
use the WiggleZ measurement of A{z = 0.6) combined with the LRG measurements of do. 2 and do. 35, indicated by bold font. 





A{z) 


dz 


Dv{z)/Dv{0.2) 






measured 


measured 


= d.o.2/dz 


= iAdzf-n- 


z = 0.2 (LRG) 


0.488 ±0.016 


0.1905 ±0.0061 




18.32 ± 0.59 


z = 0.35 (LRG) 


0.484 ± 0.016 


0.1097 ±0.0036 


1.737 ± 0.065 


10.55 ± 0.35 


z = 0.6 (WiggleZ) 


0.452 ± 0.018 


0.0692 ± 0.0033 


2.753 ±0.158 


6.65 ±0.32 



because the shape of the clustering pattern was marginahzed 
over in that analysis. 



8.1 BAO alone 

As a first step in the cosmological analysis we considered 
the constraints on cosmological parameters obtained using 
BAO data alone. Two of the distilled parameters allow the 
derivation of cosmological constraints based on only BAO 
measurements; the ratio of Dv{z) measurements and the 
acoustic parameter, A{z). 

Ratios of Dv (z) measurements are of particularly inter- 
est because they provide constraints on the cosmic expan- 
sion history using geometric information alone, independent 



of the shape of the clustering pattern and the absolute scale 
of the standard ruler. Plots of the resulting cosmological 
constraints in our two test models (ACDM and uiCDM) are 
shown in Figure 18. Using only LRG BAO data, the sin- 
gle available distance ratio Dv(0.35)/Dv(0.2) provided rel- 
atively weak constraints on cosmological parameters and, in 
particular, could not confirm the acceleration of the expan- 
sion of the universe. The addition of the second distance ra- 
tio based on WiggleZ data, Dv{0.6)/Dv{0.2), significantly 
improves these constraints. For the first time, purely geo- 
metric distance ratios from BAO measurements demonstrate 
that the cosmic expansion is accelerating: assuming the flat 
uiCDM model, a dark energy fluid with < — 1/3 is re- 
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Figure 18. Likelihood contours for cosmological parameter fits to BAO measurements using Dy{z) distance ratios. We fit two different 
models: curved cosmological-constant universes parameterized by (nm,f^A) and fl^^t dark-energy universes parameterized by {Qui,w). 
Blue contours show the constraints using the measurement of Dy (0.35) / Dy (0.2) obtained by Percival et al. (2010). Red contours display 
the new constraints in Dy (0.6) /Dy (0.2) derived using WiggleZ data. The combination of these measurements is plotted as the grey 
shaded contours. One dimensional marginalized likelihoods for fl\ and w are displayed on the right-hand side of each contour plot. In 
the flat idCDM model, the BAO distance ratios alone require accelerating expansion (w < —1/3) with a likelihood of 99.8%. Compared 
to the LRG data, WiggleZ does not favour as high values of Oa or as negative values of w. 




Figure 19. A comparison of the WiggleZ results and previous BAO measurements. The blue and green solid contours show the 
redshift z = 0.2 and z = 0.35 bins from Percival et al. (2010), using the dz parameter and including a CMB-motivated prior Qmh^ = 
0.1326 ± 0.0063. The red contour displays the fit to the WiggleZ measurement of A{z) at z = 0.6, which is independent of any prior in 
Qmh'^. We note that the degeneracy directions in these plots rotate with redshift, demonstrating the utility of combining measurements 
at different redshifts. The grey shaded region shows the combination of do. 2 and do. 35 (including the correlation). The black solid lines 
are the 1, 2 and 3-0- contours of the final result combining all these BAO measurements. The addition of WiggleZ data reduces the l-tr 
uncertainty by about 50% in Qa and about 30% in to, particularly disfavouring high values of f2A and more negative values of w. We note 
that in our final cosmological measurements we do not use a prior in Qmh^ but instead combine our results with the WMAP distance 
priors. 



quired with a likelihood of 99.8% (assuming the flat prior 
-3 < ^i; < 0). 

The improvement in the statistic comparing the 
best-fitting ACDM model [Qm, ^a) = (0.25, 1.1) to the Ein- 
stein de-Sitter model (Qn^MA) = (1.0,0.0) is Ax^ = 18, 
whilst comparing the best-fitting ACDM model to the open 



CDM model (rimMA) = (0.27,0.0) we obtain Ax^ = 8. 
Even given the extra parameter in the ACDM model, infor- 
mation criteria tests consider this level of improvement in 
X^ to be significant evidence in favour of the ACDM model 
compared to a model with no dark energy. 

Following the addition of the WiggleZ measurement, the 
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BAO data alone require accelerating cosmic expansion with 
a higher level of statistical confidence than the initial lumi- 
nosity distance measurements from supernovae that are con- 
sidered to be the first direct evidence of accelerating expan- 
sion (compare the left-hand panel of Figure 18 with Figure 6 
of Riess et al. 1998). Although these BAO measurements are 
not yet competitive with the latest supernova constraints, 
it is nevertheless reassuring that a standard-ruler measure- 
ment of the expansion of the universe, subject to an entirely 
different sot of potential systematic uncertainties, produces a 
result in agreement with the standard-candle measurement. 

The cosmological constraints from the acoustic param- 
eter A{z) arc much more constraining than the Dv{z) ratios 
because they implicitly incorporate a model for the cluster- 
ing pattern and standard ruler scale as a function of flmh^ ■ 
Figure 19 displays the resulting cosmological parameter fits 
to the WiggleZ measurement of A{z = 0.6) combined with 
the LRG measurements of c?z=o.2,o.35. For the purposes of 
this Figure we combine the measurements with a prior 
flm/i^ = 0.1326 ± 0.0063 (following Percival et al. 2010); 
in the next sub-section we use the WMAP distance priors 
instead. The improvement delivered by the WiggleZ data 
can be seen by comparing the shaded grey contour, which 
is the combination of the LRG results, with the solid black 
contours representing the total BAO constraint including 
the new WiggleZ data. The current WiggleZ dataset de- 
livers an improvement of about 50% in the measurement 
of Q,A and about 30% in the measurement of w, based on 
LSS data alone. The marginalized parameter measurements 
are = 0.25t^;g| and Oa = l.lloj (for ACDM) and 
Q.m = 0.23 ± 0.06 and w = -1.6l^;? (for wCDM). 

The BAO results continue to prefer a more negative 
dark energy cquation-of-statc or a higher cosmological con- 
stant density than the CMB or supernovae data. We explore 
this further in the following Section. 

8.2 BAO combined with CMB and SNe 

We now combine these large-scale structure (LSS) mea- 
surements with other cosmological datasets. We incorpo- 
rated the CMB data using the WMAP distance priors in 
{£a,TZ,z*) described in Section 7.1. We fitted a model pa- 
rameterized by (ilmjWjh) using flat priors 0.1 < < 0.5, 
0.5 < h < 1.0 and — 3 < w < and assuming a flat Uni- 
verse (f2k = 0). Figure 20 compares the combined LSS cos- 
mological parameter measurements to the CMB constraints 
in both the (f2m,w) and {h,w) planes, marginalizing over 
h and fim respectively. The marginalized measurements of 
each parameter arc ftm = 0.287l° o28, w = -0.982lg jg9 

and /i = o.692i;;:'(!;^^. 

We also combined these constraints with those aris- 
ing from type la supernovae data from the "Union2" com- 
pilation by Amanullah et al. (2010), which includes data 
from Hamuy et al. (1996), Riess et al. (1999, 2007), Astier 
et al. (2006), Jha et al. (2006), Wood-Vasey et al. (2007), 
Holtzman et al. (2008), Hickcn et al. (2009) and Kcssler et 
al. (2009). Using aU of these data sets (LSS-hCMB-l-SN) 
the best fitting wCDM model is {nm,w) = (0.284 ± 
0.016, -1.026 ± 0.081). This model provides a a good fit to 
the data with a minimum per degree of freedom of 0.95. 

The best-fitting parameter values based on LSS data 
alone are offset by slightly more than one standard devi- 



ation from the best-fitting parameter values of the com- 
bined LSS+CMB+SN fit. This mild tension between LSS 
and CMB-I-SN results was already evident in the Percival 
et al. (2007) measurement of BAOs in the SDSS: the LSS 
measurements prefer a slightly higher rate of acceleration, 
or a more negative equation of state of dark energy. The new 
WiggleZ results reported here serve to amplify that tension, 
although the level of discrepancy is not statistically signifi- 
cant. 

In order to investigate whether more complex models 
give a better statistical fit to these datasets we performed fits 
of two additional cosmological models. The first is a wCDM 
model in which we allow curvature to be a free parameter; 
i.e. we fit for f2m, f2k and w. The second is a «;(a)CDM 
model in which the equation-of-state of dark energy is al- 
lowed to vary linearly with scale factor (CPL parameteriza- 
tion; Chevallier & Polarski 2001, Linder 2003); i.e. we fit for 
Jim, and Wa- If these additional parameters are justified 
we should find that the minimum value of the statistic 
decreases significantly. However, we in fact find that the ad- 
dition of the extra degrees of freedom improves the value 
of the best-fitting model by at most Ax^ = 0.6 compared to 
the ACDM model, and significantly degrades the errors in 
our fitted parameters. In terms of information criteria this 
does not represent a sufficient improvement to justify the 
addition of the extra degree-of-freedom. 



9 CONCLUSIONS 

We summarize our conclusions as follows: 

• This intermediate sample of the WiggleZ Dark Energy 
Survey, containing A'^ = 132,509 galaxies, permits a convinc- 
ing detection of baryon acoustic oscillations at the highest 
redshift achieved to date, z = 0.6. We have quantified the 
baryon oscillations in both the galaxy correlation function 
and power spectrum statistics. The statistical significance of 
the correlation function detection exceeds 3-a. 

• We present the first measurement from a galaxy survey 
of the band-filtered BAO estimator wo(r) defined by Xu 
et al. (2010), which also exhibits strong evidence of BAOs, 
and we suggest some improvements in its implementation. 
The distance-scale measurements resulting from wo{r) are 
slightly less precise than those obtained by fits to the stan- 
dard correlation function, likely due to the suppression by 
the filtering function of clustering information on small and 
large scales. 

• The clustering statistics are well-fit by a non-linear 

power spectrum model including a Gaussian damping fac- 
tor which models coherent flows on 100 h^^ Mpc scales. We 
show that our results axe not sensitive to the details of the 
quasi-linear model by comparing five different implementar 
tions suggested in the literature. We use the GiggleZ N-body 
simulation to demonstrate that scalo-depondont bias effects 
in the WiggleZ galaxy distribution have an amplitude of 1% 
on 10 h^^ Mpc scales, compared to 10% for more highly- 
biased LRGs. 

• We measured the distance-redshift relation Dv{z = 
0.6) with an accuracy of about 5% using the scale di- 
lation method. The self-consistency of the measurement 
amongst the independent clustering statistics suggests that 
systematic measurement errors are not dominating these 
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Figure 20. Likelihood contours for cosmological parameter fits to LRG and WiggleZ BAO data, compared to and combined with CMB 
measurements. We have fitted for three parameters {Qm, w and h) assuming a fiat Universe. The left-hand panel of contours marginalizes 
over h and the right-hand panel of contours marginalizes over Qrn- The LSS constraints appear in blue, as dotted lines showing the Percival 
et al. (2010) results for dz with no CMB prior, and as solid lines combining that data with the WiggleZ measurement of A{z). The red 
contours show the CMB results (using the WMAP distance priors) . The grey shaded contours display the combination of BAO and CMB 
measurements. One dimensional likelihood distributions for each of the parameters are shown as the insets. The BAO constraints shown 
here are much wider than those plotted in Figure 19 because they do not assume any prior in flmh?- The improvement contributed by 
the BAO, as compared to the CMB constraints alone, can be seen by comparing the black likelihood distributions, representing the total 
constraints, with the red likelihood distributions, which are for the CMB alone. (The sharp edges of the CMB likelihood distributions 
arise due to the flat priors 0.5 < h < 1.0 and 0.1 < f2m < 0.5 we adopted). 



results. Our data best constrain the quantity A{z) = 
Dv (z) ^ ilm-ffo /'-•^i which lies perpendicular to the principal 
degeneracy direction of our contours and which is measured 
with an accuracy of 4%: we find A{z = 0.6) = 0.452 ± 0.018. 

• The distance ratios between our measurements and pre- 
vious analyses of the Luminous Red Galaxy distribution at 
redshifts z — 0.2 and z = 0.35 are consistent with a flat 
ACDM cosmological model which also provides a good fit 
to the CMB distance priors. Addition of the WiggleZ data 
allows us to estabhsh, using geometric distance ratios alone, 
that the equation-of-state of dark energy drives an accelerat- 
ing expansion (w < —1/3) with 99.8% likelihood, assuming 
a flat prior —3 < w < 0. The current WiggleZ dataset de- 
livers an improvement of about 50% in the measurement of 
f^A and 30% in the measurement of w, based on BAO data 
alone. The WiggleZ measurement confirms the mild tension 
that was previously reported between CMB and BAO mea- 



surements of the acceleration of the Universe, whereby the 
BAO data favour a slightly higher acceleration rate. 

• Cosmological parameter fits using BAO and CMB data 
are consistent with those based on current supernovae data; 
this cross-check does not yield evidence for systematic er- 
rors. Combining all current BAO, CMB and supernovae data 
we find that the best-fitting wCDM model is (Qunw) = 
(0.284 ± 0.016, -1.026 ± 0.081). These data do not justify 
the addition of another degree-of-freedom such as non-zero 
spatial curvature or an evolving equation-of-state of dark 
energy. 

When the complete BAO dataset from the WiggleZ sur- 
vey is available, we plan to split our data into redshift bins 
and explore fitting for the BAOs separately in the tangen- 
tial and radial directions, improving the cosmological con- 
straints presented in this paper. 

The WiggleZ dataset enables further tests of the cos- 
mological model, complementary to those involving baryon 
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acoustic oscillations. Two such results will be published at 
a similar time as the current paper. Firstly, we have used 
redshift-spaco distortions in the WiggleZ sample to measure 
accurately the growth rate of structure in the redshift range 
0.1 < z < 0.9, finding that this is consistent with the pre- 
dictions of the ACDM model (Blake et al. 2011a). Secondly, 
we have used the Alcock-Paczynski test to perform non- 
parametric reconstructions of the cosmic expansion history 
(Blake et al. 2011b). 
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